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ABSTRACT 

Most of the baryons in the Universe are not in the form of stars and cold gas in galaxies. 
Galactic outflows driven by supernovae/stellar winds are the leading mechanism for 
explaining this fact. The scaling relation between galaxy mass and outer rotation 
velocity (also known as the baryonic Tully-Fisher relation, BTF) has recently been 
used as evidence against this viewpoint. We use a ACDM based semi-analytic disk 
galaxy formation model to investigate these claims. In our model, galaxies with less 
efficient star formation and higher gas fractions arc more efficient at ejecting gas from 
galaxies. This somewhat counter intuitive result is due to the (observational) fact that 
galaxies with less efficient star formation and higher gas fractions tend to live in dark 
matter haloes with lower circular velocities, from which less energy is required to escape 
the potential well. In our model the intrinsic scatter in the BTF is ~ 0.15 dex, and 
mostly reflects scatter in dark halo concentration. The scatter is largely independent of 
galaxy structure because of the large radius within which galaxy rotation velocities are 
measured. The observed scatter, equal to ~ 0.24 dex, is dominated by measurement 
errors. The best estimate for the intrinsic scatter is that it is less than 0.15 dex, and 
thus our ACDM based model (which does not include all possible sources of scatter) 
is only just consistent with this. Future observations of the BTF scatter could be 
made with a more stringent measurement of the intrinsic scatter, and thus provide 
a strong constraint to galaxy formation models. In our model, gas rich galaxies, at 
fixed virial velocity (V™), with lower stellar masses have lower baryonic masses. This 
is consistent with the expectation that galaxies with lower stellar masses have had 
less energy available to drive an outflow. However, when the outer rotation velocity 
(Vflat) is used the correlation has the opposite sign, with a slope in agreement with 
observations. This is due to the fact there is scatter in the relation between Va a t and 
Vvir- In summary, contrary to some previous claims, we show that basic features of the 
BTF are consistent with a ACDM based model in which the low efficiency of galaxy 
formation is determined by galactic outflows. 

Key words: galaxies: formation - galaxies: fundamental parameters - galaxies: haloes 
- galaxies: spiral - galaxies: kinematics and dynamics 



1 INTRODUCTION 

The cosmic baryon fraction is extremely well determined 
from observations of the CMB plus other cosmological 
probes, with the latest results from WMAP finding /b ar = 
n b /fi m = 0.167 ±0.004 (Komatsu et al. 2011). However, on 
galaxy scales a significant fraction of the baryons are "miss- 
ing". Stars and cold gas in galaxies account for just ~ 8% of 
the cosmic baryons (e.g., Bell et al. 2003; Fukugita & Pee- 
bles 2004; Read & Trentham 2005), while hot intracluster 
gas accounts for just ~ 4% (Fukugita & Peebles 2004). 

The vast majority of the cosmic baryons (~ 88%) are 
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thought to be in the form of hot gas in the haloes of 
galaxies or between galaxies in the so called warm-hot- 
intergalactic medium (WHIM) at temperatures between 10 5 
and 10 7 K (Cen & Ostriker 1999). However, only a fraction 
of these baryons have been detected (e.g., Bregman 2007, 
Shull et al. 2011), and the amount of baryons that reside in 
hot haloes around the Milky Way and other nearby galax- 
ies is a subject of current debate (e.g., Grcevich & Putman 
2009; Anderson & Bregman 2010). 

This raises the question: Why are most of the cosmic 
baryons in hot haloes or the WHIM? There are two basic 
answers: 1) Most of the baryons accreted into galaxies and 
were then expelled (into haloes or the WHIM) by feedback 



© 2012 RAS 



2 Button 



processes (stellar and/or AGN); or 2) Most baryons never 
accreted into galaxies in the first place. 

The answer to this question is of interest beyond the 
realm of baryon accounting. Outflows have been invoked 
to explain a number of apparent problems with galaxy for- 
mation in a ACDM context. These include the predicted 
central density cusps, which are not observed, but can be 
softened with galactic outflows (e.g., Read & Gilmore 2005; 
Mashchenko et al. 2006; Governato et al. 2010; Maccio 
et al. 2012), and the excess of low-angular momentum mate- 
rial which needs to be removed in order to produce bulgeless 
disk galaxies with exponential density profiles (e.g., Mailer 
& Dekel 2002; Dutton 2009; Governato et al. 2010; Brook 
et al. 2011). Thus if galactic outflows are not the explanation 
of why galaxy formation is so inefficient, other mechanisms 
will need to be found to reconcile ACDM with observations 
on galaxy scales. 

A clue to the origin of the missing baryons comes from 
the fact that the galaxy formation efficiencjQ is not a con- 
stant. It is observed to peak at egf ~ 30% in haloes of 
mass M v i r ~ 1O 12 M0, and declines to both higher and lower 
masses (e.g., Hoekstra et al. 2005; Conroy & Wechsler 2009; 
Dutton et al. 2010; Moster et al. 2010; More et al. 2011). In 
high mass haloes M v i r > 1O 12 M0, cooling is progressively 
more inefficient (e.g., Blumenthal et al. 1984), which results 
in the correct qualitative trend. However, additional heat- 
ing mechanisms (such as AGN feedback) are needed in order 
to reproduce the rapid decline in galaxy formation efficiency 
with increasing halo mass (e.g., Croton et al. 2006). In haloes 
of mass 10 10 < M v i r < 10 12 Mq most of the cosmically avail- 
able baryons should accrete onto central galaxy (e.g., Blu- 
menthal et al. 1984; Keres et al. 2009). In this mass range, 
galactic outflows driven by supernovae or stellar feedback 
are the leading mechanism for explaining why the galaxy 
formation efficiencies are so low, and decline with decreas- 
ing halo mass. 

Galactic outflows appear ubiquitous in galaxies that are 
undergoing, or have recently undergone, significant star for- 
mation (e.g., Shapley et al. 2003; Martin 2005; Tremonti 
et al. 2007; Weiner et al. 2009; Rubin et al. 2010; Steidel 
et al. 2010). At least some of the outflowing gas is observed 
to be moving faster than the escape velocity of the halo. 
However, measuring outflow mass rates is challenging, and 
at present it is not clear how much mass is actually re- 
moved (e.g., Rubin et al. 2010). Thus, while galactic out- 
flows undoubtedly exist, their role in determining the bary- 
onic masses of galaxies is unclear. 

The scaling relations between baryonic mass, outer ro- 
tation velocity and gas fraction have been used as arguments 
against galactic outflows being the explanation for the ob- 
served low galaxy formation efficiencies. Anderson & Breg- 
man (2010) argue that outflows should result in a negative 
correlation between galaxy mass and stellar mass at fixed ve- 
locity, while no such correlation is observed. They cite this 
as strong evidence against galatic outflows. McGaugh (2012) 
shows that the efficiency of outflows needs to be higher in 
galaxies with higher gas fractions and lower past average star 



1 We define the galaxy formation efficiency as the fraction of the 
cosmically available baryons, /b ar > that end up as stars and cold 
gas in a galaxy: e GF = M ga i/(/b a rMvir)- 



formation rates. McGaugh (2012) argues that in the context 
of feedback models this is apparently puzzling because of the 
notion that galaxies with more star formation should have 
more energy to drive an outflow. Hence the galaxies that are 
most efficient at removing baryons are expected to have the 
highest star formation efficiencies and lowest gas fractions. 
As we show below the resolution of this puzzle is the fact 
that it is not just the amount of star formation that deter- 
mines the efficiency of feedback. The depth of the potential 
well is also important — it is much easier to remove baryons 
from lower mass haloes. In addition McGaugh (2011, 2012) 
has argued that the scatter in the BTF is consistent with 
being zero, which is hard to explain in a ACDM context. 

It should be noted that until the details of star for- 
mation and feedback are understood it will not be possible 
to talk of definitive predictions for the BTF in the ACDM 
paradigm. The question we can ask at the present time is 
whether the properties of the BTF can be reproduced in 
a ACDM context using plausible models for star formation 
and feedback. In this paper we address this question using 
the semi-analytic disk galaxy formation model of Dutton & 
van den Bosch (2009). This paper is organized as follows: 
In §2 we give a brief outline of the galaxy formation model; 
In §3 we discuss the correlations between ejection efficiency, 
galaxy velocity and gas fraction; In §4 we discuss the scatter 
in the baryonic Tully-Fisher relation; A summary is given 
in §5. 



2 GALAXY FORMATION MODEL 

In this paper we use a sample of galaxies generated from a 
ACDM based semi-analytic galaxy formation model (Dut- 
ton & van den Bosch 2009). Very briefly, this model fol- 
lows the evolution of resolved disks of gas and stars in- 
side smoothly growing dark matter haloes. The mass accre- 
tion histories, dark halo structure and angular momentum 
are consistent with cosmological N-body simulations (Bul- 
lock et al. 2001a,b; Wechsler et al. 2002; Sharma & Stein- 
metz 2005; Maccio et al. 2007, 2008). We create a Monte- 
Carlo sample of 2000 galaxies by uniformly sampling halo 
masses from \og(AI v i T /h~ 1 Mq) = 9.5 to 13.5, and adding 
log-normal scatter in halo concentrations (<Ti Ilc = 0.25), spin 
parameters (<j\ n \ = 0.50), and angular momentum profile 
shapes (<ri na = 0.25). The feedback efficiency (efb = 0.5) 
and angular momentum losses are tuned to match the galaxy 
formation efficiency vs halo mass and galaxy specific angular 
momentum vs halo mass relations (see Dutton & van den 
Bosch 2012). 

The main limitations of this model are discussed in Dut- 
ton & van den Bosch (2009). Here we highlight two effects 
that are most likely to impact the BTF - the assumptions 
of smooth mass accretion and that outflow gas does not re- 
turn to the galaxy. While smooth accretion is expected to 
dominates the build up of spiral galaxies in ACDM cosmolo- 
gies, deviations from this in the form of minor and major 
mergers are an unavoidable feature. In our outflow model 
we only consider winds that escape the halo, and do not re- 
turn. It is likely that some of the ejected gas will return to 
the halo at a later time (e.g., Oppenheimer et al. 2010). In 
addition, it is possible for gas to escape the disk, but not 
the halo, and also be reaccreted at a later time - which is 
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Figure 1. Relations between mass loading factor, outer galaxy circular velocity and gas-to-stellar mass ratio for a ACDM based galaxy 
formation model with galactic outflows. Left: Feedback is more efficient at ejecting baryons from lower mass galaxies, despite their lower 
star formation efficiencies. Middle: Feedback is more efficient at ejecting baryons from galaxies with higher gas-to-stellar mass ratios. 
Right: Galaxies with lower velocities have higher gas fractions, both in models (points, solid line) and observations from Stark et al. (2009) 
and McGaugh (20f2) (open points). This latter correlation explains why middle and left panels are mirror images of each other. 



known as a galactic fountain (e.g., Brook et al. 2012a). All 
of these effects create scatter in the mass accretion histories 
of baryons and dark matter onto galaxies, which one would 
nominally expect to result in more scatter in the BTF. 



3 EJECTION EFFICIENCY 

We define the ejection fraction as the fraction of the cos- 
mically available baryons that have been ejected from the 
galaxy: 

^cjcct 

= M e ject/(/barM v j r ). Similarly we define the 
stellar mass fraction as the fraction of the cosmically avail- 
able baryons that end up in stars (or stellar remnants): 
tstar = A/star/(/barM v ir). Note that due to return of gas 
from stars into the ISM, e sta r is less than the integral of the 
star formation rate for any given galaxy. The ratio of these 
two quantities is the mass loading factor: rj — e c ject/e s tar, 
which is a way to parametrize how "efficient" feedback is 
in ejecting gas from a galaxy with a given amount of star 
formation. 

The left panel of Fig. [1] shows the mass loading fac- 
tor vs outer galaxy circular velocity, VflaJE for our model. 
Galaxies with lower circular velocities have high mass load- 
ing factors (i.e., lower mass galaxies are more efficient at 
ejecting their baryons). As in the observations (McGaugh 
2012) the effect is not subtle: galaxies with Va a t ~ 40kms _1 
have T) ~ 100, while galaxies with Va at ~ 200 km s have 
r\ ~ 1. The scaling between mass loading factor and circular 
velocity 77 oc V~ 2 (dashed line) results from our assumptions 
of energy conservation and that the outflow moves at the lo- 
cal escape velocity. We note that momentum driven winds 
(e.g., Murray et al. 2005) are expected to result in 77 oc V~ x , 
while constant wind velocity models (e.g., Springel & Hern- 
quist 2003) result in a constant mass loading factor. Thus 
as long as the wind velocity scales with the galaxy velocity, 
the lower energy/momentum input from stars/supernovae 



2 In our model we measure Vfl a t at the radius at a radius which 
encloses 80% of the cold gas. 



in lower mass galaxies is more than compensated by the 
shallower potential wells. 

The middle panel of Fig.[T]shows the mass loading factor 
vs gas-to-stellar mass ratio, JV/g as /M star . Where M gas is the 
mass in cold gas (atomic and molecular). This shows that 
galaxies with higher gas fractions have higher mass load- 
ing factors, and are thus more "efficient" at ejecting their 
baryons, in qualitative agreement with observations (Mc- 
Gaugh 2012). A comparison between the left and middle 
panels shows that there is less scatter in the relation be- 
tween mass loading factor and galaxy velocity, than between 
mass loading factor and gas fraction. This suggests the rela- 
tion between mass loading and galaxy velocity (left panel) is 
more fundamental than the relation between mass loading 
and gas fraction (middle panel). The reason galaxies with 
higher gas fractions have higher mass loading is simply a re- 
sult of the anti-correlation between gas fraction and galaxy 
velocity (right panel): (M gaa /M s t ar ) oc Vg"*' 6 . As before, the 
solid line shows the median of the model (grey points) which 
is in good agreement with the observations (open circles) 
from Stark et al. (2009) and McGaugh (2012). 

This raises the question: Why do lower mass galaxies 
have, on average, higher gas fractions? Observationally we 
know that lower mass galaxies are on average less dense 
(e.g., Kauffmann et al. 2003), and that lower density galaxies 
are less efficient at turning gas into stars (Kennicutt 1998). 
This leads naturally to higher gas fractions in lower mass 
galaxies. However, to reproduce the observations in detail 
requires, in addition to the standard Schmidt-Kennicutt star 
formation law, a threshold density for star formation (van 
den Bosch 2000). On the theory side, a correlation between 
galaxy density and galaxy mass occurs naturally in a ACDM 
context. The simplest disk galaxy formation model (in which 
the galaxy formation efficiency and spin parameters are con- 
stant) results in disk sizes Rd oc M^ 3 and thus disk densities 

E d oc M g ^ 3 (e.g., Mo, Mao, & White 1998). Including out- 
flows typically results in shallower size- mass relations (e.g., 
Dutton & van den Bosch 2009), and thus an even stronger 
mass - gas density relation. 
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Figure 2. Relations between baryonic (Mg a i) and stellar (M s tar) mass with virial (V v ir) an d outer galaxy (Vfl a t) velocity for a ACDM 
based galaxy formation model (grey points, solid line) and observations from Stark et al. (2009) and McGaugh (2012). The slopes of 
all three relations deviate from the /bar-^vir °c KiL scaling for ACDM haloes, which is indicated by the dotted lines. In the models the 
scatter increases from 0.12 dex for the M ga j vs V v i r relation to 0.15 dex for M ga i vs Vfl a t (BTF) and 0.22 dex for M B t&r vs Vfl a t (STF). 



4 THE BARYONIC TULLY-FISHER 
RELATION 

The baryonic Tully-Fisher relation (BTF) is the relation be- 
tween the baryonic mass of a galaxy, M ga i (stars and cold 
gas), and the rotation velocity at large galactic radii, typ- 
ically referred to as Vflat- It is an extension of the original 
Tully-Fisher relation which is a correlation between HI line- 
width and galaxy luminosity (Tully & Fisher 1977). The 
BTF was first studied by McGaugh et al. (2000), and sub- 
sequently by numerous authors, both observationally (e.g., 
Bell & de Jong 2001; McGaugh 2005; Geha et al. 2006; Avila- 
Reese et al. 2008; Begum et al. 2008; Stark et al. 2009; Tra- 
chternach et al. 2009; Gurovich et al. 2010; Hall et al. 2011; 
Catinella et al. 2012); and theoretically (e.g., Dutton & van 
den Bosch 2009; de Rossi et al. 2010; Dutton et al. 2011; Pi- 
ontek & Steinmetz 2011; Trujillo- Gomez et al. 2011; Brook 
et al. 2012b). In this section we discuss aspects of the slope 
and scatter in the context of ACDM. 

It has been argued that the "true" BTF should include 
the contribution of ionized gas (Gnedin 2011). However, this 
is both hard to measure and hard to model. In this paper 
we restrict the baryonic mass (of both our models and the 
data) to be that of the stars and cold gas in a galaxy. Since 
the vast majority of the ionized gas will be at radii beyond 
the HI disk, the baryonic mass that we use, M ga i, is close 
to the baryonic mass within the HI radius of the galaxy. As 
such, Mgai, is expected to be more strongly correlated to 
Vflat than the total baryonic mass inside the virial radius 
of the dark matter halo, Mb ar - Indeed, in our model the 
relation between Afbar and Vflat has a slightly larger scatter 
than the regular BTF. Likewise the relation between Mbar 
and Vvir has smaller scatter than the relation between Mb ar 
and Vflat. 



4.1 Slope and zero point of the BTF in LCDM 

In ACDM the underlying origin of the BTF is the M vir oc 
Vvi T relation of dark matter haloes. This relation has no scat- 
ter by definition. Predicting the BTF from this relation re- 
quires understanding how baryonic mass is related to virial 



mass (f 9 = Mgai /M v ir); and how galaxy velocity is related 
to virial velocity: fv = Vkat/Vdr. Galaxy formation effi- 
ciencies depend on the details of gas cooling, feedback and 
recycling - none of which can be predicted from first prin- 
ciples. The relation between galaxy and halo velocities is 
better constrained (thanks to cosmological N-body simula- 
tions), but it too depends on a couple of unknown factors: 
the response of the halo to galaxy formation, and the galaxy 
formation efficiency. Because the enclosed dark matter frac- 
tions increase with increasing galactio-centric distance, these 
two unknowns are minimized by using velocities measured 
at large radii. Given there are no unique predictions for the 
BTF in ACDM cosmologies, and are unlikely to be so for 
some time, its main utilization is likely to be as a constraint 
to galaxy formation models, and in particular models for 
feedback. Indeed the feedback efficiency and angular momen- 
tum loss in our model have been tuned to match the galaxy 
formation efficiency as a function of halo mass. 

4.2 Scatter in the BTF 

The intrinsic scatter in the BTF is observed to be small 
and is consistent with zero. This is a potential problem for 
ACDM (McGaugh 2011, 2012). Fig. shows three Tully- 
Fisher relations from our galaxy formation model: Baryonic 
mass vs virial velocity (left); Baryonic mass vs galaxy veloc- 
ity (BTF, middle); and Stellar mass vs galaxy velocity (STF, 
right). All relations have small, but non-negligible, scatter of 
~ 0.1 — 0.2 dex, with the BTF having a scatter of 0.15 dex. 
There is also a velocity dependence to the scatter in all three 
relations. For the BTF the scatter ranges from ~ 0.20 dex 
at Vflat ~ 30kms _1 to ~ 0.11 dex for V fla t ~ 200 km s" 1 . 

What is the source of the scatter? In our models the 
majority (~ 73%) of the scatter in the BTF comes from 
variation in the concentration of the dark matter halo, 
which mostly effects Vrat/Vvir- The scatter in halo con- 
centrations is constrained by cosmological simulations (e.g., 
Maccio et al. 2008), and so there is little freedom to change 
this in the context of ACDM. The remainder of the scatter 
comes from variation in the halo spin parameter, which in 
our model effects the galaxy formation efficiency. We note 
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Figure 3. Correlations between residuals of the baryonic TF and stellar mass TF relations for gas rich galaxies (M gas /M s t a r > 1)- At 
fixed dark matter halo virial velocity, V v i T , lower stellar masses result in higher baryonic masses (left panel), as might be expected by 
the lower energy available to drive an outflow. However, when computed at fixed galaxy rotation velocity, Vfl at , this correlation has the 
opposite sign (middle panel), with a slope in agreement with observations (right panel). 



that while variation in the halo spin parameter is the pri- 
mary source of variation in the distribution of baryons in 
galaxies (i.e., galaxy sizes), this does not significantly effect 
Vflat because we are measuring circular velocity at a radius 
which encloses most of the baryons. Since there are likely 
other sources of scatter that are not taken into account in 
our models, we expect that our models provide a lower limit 
to the BTF scatter in ACDM. 

What ts the intrinsic scatter of the observed BTF? 
The observed scatter in the BTF is ~ 0.24 dex (McGaugh 
et al. 2011; Hall et al. 2011), and is dominated by measure- 
ment errors on baryonic masses (McGaugh 2011; Foreman & 
Scott 2011). Unfortunately, the errors have uncertainties, so 
we currently do not have a robust measurement of the intrin- 
sic scatter. Nevertheless, McGaugh (2012) finds the intrinsic 
scatter to be < 0.15 dex, which our model is just consistent 
with. 

Thus the intrinsic scatter of the BTF has the poten- 
tial to be a powerful constraint on galaxy formation models. 
In order to make progress a BTF sample needs to be con- 
structed for which the measurement errors are controlled 
to be smaller than the scatter one is trying to measure. 
There are two primary sources of measurement errors: dis- 
tance uncertainties and stellar mass uncertainties. Rotation 
velocity errors are typically small when resolved HI rota- 
tion curves are obtained, and gas masses can be measured 
reliably. Distance uncertainties can be minimized by using 
galaxies at large enough distances such that peculiar veloc- 
ities are not important (i.e., > 100 Mpc). But obtaining 
resolved HI rotation curves for such galaxies is currently a 
challenge, and may have to wait until the next generation 
radio telescopes. HI line widths or Ha rotation curves can 
be measured for large (~ 1000's) samples of galaxies (e.g., 
Courteau et al. 2007; Hall et al. 2012), but these are not 
as straightforward to interpret as Vflat ■ For gas poor spiral 
galaxies, stellar mass uncertainties are likely to remain the 
largest source of error, but it is possible for them to be ac- 
curate to ~ 0.05 dex (Gallazzi & Bell 2009), assuming one 
knows the form of the stellar initial mass function. For gas 
rich galaxies, the contribution of molecular gas is the largest 
source of error, but these galaxies are thought to be domi- 



nated by atomic gas which can be reliably measured. Thus, 
in principle, it should be possible to control measurement 
uncertainties to less than 0.1 dex, and thus determine if the 
intrinsic scatter of the BTF is smaller than the 0.15 dex of 
our model. 



4.3 Correlation between BTF and STF 

Anderson & Bregman (2010) used the scatter in the BTF 
to test the idea that galactic outflows are responsible for 
the low observed galaxy formation efficiencies. The idea is 
that for a galaxy in a given dark matter halo, a galaxy that 
forms more stars will have more energy /momentum avail- 
able to drive an outflow, and thus should remove a higher 
fraction of its baryons. Thus at fixed velocity, there should 
be a negative correlation between the baryonic mass and 
stellar mass. Anderson & Bregman (2010) found no such cor- 
relation, and thus argued this was strong evidence against 
galactic outflows being responsible for the low galaxy for- 
mation efficiencies. 

We investigate the validity of this reasoning using BTF 
and STF relations from a semi-analytic model. Fig. shows 
the correlations between the mass residuals of these BTF 
and STF relations. We show results using TF relations con- 
structed using both virial velocity of the dark matter halo, 
Vvir, (which is not observable for individual galaxies), and 
the velocity in the outer part of a galaxy, Vfl at , (which is 
observable for individual galaxies) . We split the models into 
gas rich and gas poor, since gas poor galaxies are (trivially) 
expected to have a positive correlation between baryonic 
mass and stellar mass (Fig. [4}. The upper left panel shows 
that at fixed V v i r , gas rich galaxies do indeed have a neg- 
ative correlation between baryonic mass and stellar mass. 
However, at fixed Vfl at , the correlation has the opposite sign 
(middle panel), with a slope in good agreement with obser- 
vations from Stark et al. (2009) and McGaugh et al. (2012) 
(right panel). This change in slope is a result of scatter in 
the relation between Vfl a t and V v ir- Thus the simple reason- 
ing used by Anderson & Bregman (2010) to argue against 
outflow models, while correct in principle, is not valid in 
practice. 
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Figure 4. As Fig.[3]but for gas poor galaxies (M gas /M s tar < 1)- There is a strong correlation between residuals of the baryonic TF and 
stellar TF relations, in both models and data. This is a trivial consequence of the fact that M s tar — M~~\ in these galaxies. 



5 SUMMARY 

We have used a ACDM based galaxy formation model to 
investigate the observable signatures of galactic outflows in 
the baryonic Tully-Fisher relation (BTF). We summarize 
our results as follows: 

• Observations indicate that galaxies with lower star for- 
mation efficiencies and higher gas fractions have higher ejec- 
tion efficiencies (e.g., McGaugh 20ff). We show that these 
trends can be explained by energy driven feedback models. 

• In our model feedback is more efficient in galaxies with 
lower circular velocities (and shallower potential wells). This 
is in spite of the significantly lower star formation efficiencies 
in lower velocity galaxies. 

• In our model, as well as observations, lower velocity 
galaxies have higher gas fractions. This results in a (non- 
causal) correlation between ejection efficiency and gas frac- 
tion, such that ejection efficiencies are higher in galaxies 
with higher gas fractions. 

• Lower mass galaxies are predicted to have higher gas 
fractions. This is a generic prediction for galaxy formation 
in ACDM, which is strengthened by feedback. 

• The scatter our model BTF is ~ 0T5 dex, and is mostly 
due to variations in the dark matter concentration parame- 
ter. While this scatter is significantly smaller than the ob- 
served scatter of ~ 0.24 dex, most of the observed scatter 
is due to measurement uncertainties. McGaugh (2012) finds 
the intrinsic scatter in the BTF is < 0T5 dex, which our 
model is only just consistent with. 

• In principle, future observations of the BTF could be 
made where the measurement uncertainties are controlled to 
less than 0.1 dex, and thus to provide stringent constraints 
to the intrinsic scatter, and ACDM based galaxy formation 
models. 

• In our model, at fixed virial velocity (Vvir), gas rich 
galaxies with lower stellar masses have higher baryonic 
masses. This is consistent with the idea that less star for- 
mation should result in less energy (or momentum) available 
to drive an outflow (e.g., Anderson & Bregman 2010). How- 
ever, at fixed galaxy velocity (Vflat), the correlation is of the 
opposite sign due to the scatter in Vaat/VW. Furthermore 
the slope of the correlation in our model is in agreement 
with observations. 



In summary, we find that there is currently no conflict 
between the observed baryonic Tully-Fisher relation (slope, 
scatter, and residual correlations) and predictions of ACDM 
based models in which galaxy formation efficiencies are de- 
termined by galactic outflows. 
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